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Abstract 


It is shown that a spin system with long range interactions can be converted 
into a chaotic dynamical system that is differentiable and low-dimensional. The 
thermodynamic limit of the spin system is then equivalent to studying the long 
term behavior of the dynamical system. Cycle expansions of chaotic systems 
(expansion of the Fredholm determinant) are then used to study the thermo¬ 
dynamic limit. By considering the smooth dynamical system, it is possible to 
converge to the thermodynamic limit faster than with transfer matrices. 
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Cycle expansions are an efficient method to compute the properties of er- 
godic systems [1]. They have been successfully applied to chaotic dynamical 
systems [1, 2], spin systems in statistical mechanics [3, 4], and to the semiclassi- 
cal quantization of chaotic systems [5, 6]. In a cycle expansion, periodic orbits of 
the system are defined and each orbit is given a weight (such as the stability of 
an orbit, or the energy of a configuration). The weights are then combined into 
an expansion that uses shorter orbits to estimate the weight of longer orbits. 
The more orbits added to the expansion, the more accurate the estimates. The 
great advantage of cycle expansions is that analytic understanding of the system 
(on symmetries, or on partly integrable regions, for example) can be combined 
with numerical results into a rapidly converging expansion. In the case of dy¬ 
namical systems what is expanded is the Fredholm determinant of the operator 
that determines the invariant measure of the system — the Perron-Frobenius 
operator. This operator is defined in such a way that the determinant has no 
singularities, which leads to expansions that converge faster than any exponen¬ 
tial in the length of the longest periodic orbit. For statistical mechanics the 
transfer matrix is very similar to the Perron-Frobenius operator, but it leads to 
a determinant with singularities. And in quantum mechanics it is not known, 
in general, how to construct the operator equivalent to the Perron-Frobenius or 
transfer matrix (although Cvitanovic and Vattay have solved the semiclassical 
version of this problem [7]). 

I would like to show that it is possible to construct a Fredholm determi¬ 
nant for spin systems. I will do this by transforming the spin system into a 
low-dimensional dynamical system and then constructing its Fredholm deter¬ 
minant. Using the known results of chaotic systems, this means that expan¬ 
sions for thermodynamic quantities will converge faster (both qualitatively and 
quantitatively) than the exponentially fast convergence obtained using trans¬ 
fer matrices. There are many advantages of carrying out this procedure for a 
spin system: it allows the use of differentiable methods, before difficult due to 
the discrete nature of spins; and it gives a general framework in which to con¬ 
struct other Fredholm determinants, as there are few restrictions on the spin 
interaction. 

The use of methods from statistical mechanics of spin systems in the study of 
dynamical systems — the thermodynamic formalism — is now well established. 
It was pioneered by Sinai [8], Bowen [9], and Ruelle [10] in the early seventies 
and it had its great success in physics after its application to the Feigenbaum 
period doubling route to chaos [11, 12]. In contrast, the methods of dynamical 
systems have not been used to study spin systems. With the geometrization 
of spin systems this reverse route becomes possible. As a first application, 
geometrization provides a counter-example to the belief that, aside from exactly 
solved models, the thermodynamic limit can only be approached exponentially 
fast. 
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The basic idea in the geometrization of spin systems is not to consider the 
shift of spin configurations as the dynamics, but instead to consider the dynam¬ 
ics of a few observables (such as the interaction energy or the magnetization). 
To explain how this is done, I will review the transfer operator (a generalization 
of the transfer matrix) which acts on observables (functions of configurations). 
Then the £ operator will be introduced, similar to the transfer operator, but 
acting in a more restricted space, that of functions that depend only on the 
effective field at a point. This difference turns out to be crucial. To study 
the spectrum of the £ operator one uses the techniques of cycle expansions. 
A difficulty in computing the trace of is solved by considering a functional 
derivative involving the interaction. As an example of the advantages of using 
the £ operator, I show that the thermodynamic limit is approached faster than 
any exponential. 

To geometrize spin systems, the interactions are assumed to be translation- 
ally invariant and smooth. The spins will only assume a finite number of 
values. For simplicity I will take the interaction cf) among the spins to depend 
only on pairwise interactions, 

4)(p) =4 )(Po,P 1,P2,...) = JoPO + ^5cro,cr,Jl(^) (1) 

Tl>0 

and the possible values of a spin Pti to be in CIq = {-f, —For the one¬ 
dimensional Ising model, Jo is the external magnetic field and Ji (n) = 1 if n = 1 
and 0 otherwise. For an exponentially decaying interaction Ji(rL) = 

The transfer operator T was introduce by Kramers and Wannier [13] to 
study the Ising model on a strip and concocted so that the trace of its n-th 
power is the partition function Zn of system when one of its dimensions is n. 
The method can be generalized to deal with any finite-range interaction. If the 
range of the interaction is L, then T is a matrix of size 2^ x 2^. The longer 
the range, the larger the matrix. When the range of the interaction is infinite 
the matrix definition becomes cumbersome, and it is simpler to define the T 
operator by its action on an observable g. Just as the observables in quantum 
mechanics, g is a function that associates a number to every state (configuration 
of spins). The energy density and the average magnetization are examples of 
observables. From this equivalent definition one can recover the usual transfer 
matrix by making all quantities finite range. For a semi-infinite configuration 

p = {Pfl, Pi,.. 

Tg(p) =g(+Vp)e-P*'+'^^> + g(-Vp)e-P*'-'^^> . (2) 

By -b V p I mean the configuration obtained by prepending -b to the beginning 
of p resulting in the configuration {-b, Pq, Pi ,..When the range becomes 
infinite, the trT’^ is infinite and there is no longer a connection between the 
trace and the partition function for a system of size n (this is a case where 
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matrices give the wrong intuition). Ruelle [14] generalized the Perron-Frobenius 
theorem and showed that even in the case of infinite range interactions the 
largest eigenvalue of the T operator is related to the free-energy of the spin 
system and the corresponding eigenvector is related to the Gibbs state. By 
applying T to the constant observable u that returns 1 for any configuration, 
the free energy f is computed as 

— |3f(|3) = lim -^InllT’^ull . (3) 

TH(X> TL 

To construct a smooth dynamical system that reproduces the properties 
of T, one uses the phase space reconstruction technique of Packard et al. 
[15] and Takens [16], and introduces a vector of state observables x(c 7 ) = 
{x] (p),... ,xd(c 5 ')}. To avoid complicated notation I will limit the discussion 
to the example x(p) = {x+(p),x_(p)}, with x+(p) = 4)(-l- V p) and x_(p) = 
4)(— V p). We also restrict all observables to be functions only of the state 
observables {xic(p)}. That is, the only accepted observables are those g for 
which, for all configurations p, there exist an analytic function G such that 
G(xi (p),... ,xd(p)) = g(p). This at first seems a severe restriction as it may 
exclude the eigenvector corresponding to the Gibbs state. It can be checked 
that this is not the case by using the formula given by Ruelle [17] for this 
eigenvector. A simple example where this formalism can be carried out is for 
the interaction 4 )(p) with pairwise exponentially decaying potential Ji(rL) = 
(with |q| < 1). In this case 4 )(p) = ^ti>o ^a-o.o-,, and the state observables 
are x+(p) = ^ti>o ^+.crnand x_(p) = ^ti>o These two observables 

always add up to the same constant. In this case the observable x+ gives the 
energy of - 1 - spin at the origin, and x_ the energy of a — spin. 

As an example, assume that the phase space of the dynamical system is 
given by the two observables {x+ (p), x_(p)} = V p), 4)(— V p)}. In terms 
of these observables the transfer operator can be re-expressed 

TG(x(p))= ^ G (x+(ri V p) ,x_ (ri V p)) ^ ^ 4 ^ 

T|=+.- 

In this equation the only reference to the configuration p is when computing the 
new observable values x+(ri V p) and x_(ri V p). The iteration of the function 
that gives these values in terms of x+(p) and x_(p) is the dynamical system that 
will reproduce the properties of the spin system. For the simple exponentially 
decaying potential this is given by two maps, F+ and F_, that give the new 
value of the observables in terms of the old. The map F+ takes {x+(p),x+(p)} 
into {x+(-F V p),x_(-F V p)} which is {q(1 -F x+),qx_} and the map F_ takes 
{x+,x_} into {qx+, q(1 -Fx_)}. In a more general case we have maps Ft| that take 
x(p) I—> x(ri V p), and a new operator £ 

£G(x)‘iVG(x(p))= Y_ G(F,(x))e-P"" . (5) 

T| = [+.-l 
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Notice that all dependencies on a have disappeared — if we know the value of 
the state observables x the action of £ on G can be computed. 

The basic idea for the dynamics will be to choose the maps Ft| so that one 
of the state variables is the interaction energy. One can consider the two maps 
F+ and F_ as the inverse branches of a map f, that is, f^’(x) = {T^+(x),F_(x)}. 
Studying the thermodynamics of the interaction cf) is equivalent to studying the 
long term behavior of the orbits of the map f, achieving the transformation of 
the spin system into a dynamical system. 

Unlike the original transfer operator, the £ operator — acting in the space 
of observables that depend only on the state variables — is of trace class (its 
trace is finite). The finite trace gives us a chance to relate the trace of 
to the partition function of system of size n. We can do better. As most 
properties of interest (thermodynamics, fall-off of correlations) are determined 
directly from its spectrum, we can study instead the zeros of the Fredholm 
determinant det(1 — z£) by the technique of cycle expansions developed for 
dynamical systems [2]. A cycle expansion consists of finding a power series 
expansion for the determinant by writing det(1 — z£) = exp(trln(1 — z£)). The 
logarithm is expanded into a power series and one is left with terms of the form 
tr£’^ to evaluate. For evaluating the trace, the £ operator is equivalent to 


£G(x) 


dy5(y-f(x))e-P«G(y) 

R 


from which the trace can be computed: 


tr£’" 


L 

x=f°" (x) 


g-|3H(x) 

det (1 — 0xf°'^(x)) 


(6) 

( 7 ) 


with the sum running over all the fixed points of f®’^ (all spin configurations of a 
given length). Here f®’^ is f composed with itself n times, and H(x) is the energy 
of the configuration associated with the point x. In practice the map f is never 
constructed and the energies are obtained directly from the spin configurations. 
The problem has now been reduced to a dynamical system and all the techniques 
applicable to the study of Fredholm determinants of dynamical systems can be 
used. 

To compute the value of tr£’^ we must compute the value of 9xf®’^; this 
involves a functional derivative. To any degree of accuracy a number x in the 
range of possible interaction energies can be represented by a finite string of 
spins e as X = 4)(+, eo, —,...). By choosing the sequence e to have a 

large sequence of spins —, the number x can be made as small as needed, so in 
particular we can represent a small variation 5 = 4)(ri). As x+(e) = 4)(-F V e), 
it can be computed by applying the map f to 4)(e). From the definition of a 
derivative we have: 


9xf(x) 


lim 

m—)oo 


4)(e Vri*™*) - 
4)(ri(™> 


4)(e) 


(8) 
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where t]*™* is a sequence of spin strings that make 5 smaller and smaller. By 
substituting the definition of cf) in terms of its pairwise interaction J(rL) = 
and taking the limit for the sequences t]*™* = {+, ,riTR+i, T^m+ 2 , ■ ■ •} one 

obtains the limit is a if y = 1, 1 if y < 1, and 0 if y > 1. It does not depend on 
the positive value of s. 

The manipulations have up to now assumed that the map f is smooth. If the 
dimension D of the embedding space is too small, f may not be smooth (D = 
oo cannot be ruled out). Determining under which conditions the embedding 
is smooth is a complicated question related to the problem of embedding an 
ergodic systems into classical systems [18]. There is no known complete solution 
to the problem, but in the case of spin systems with pairwise interaction it is 
possible to give a rule. If the interaction is of the form 

Ti>1 k 

where pi^ are polynomials and |qic| < 1, then for y > 1 one can choose a 
dynamical system with dimension D equal to the sum of all degrees, that is 
Y_ deg(pic). The state observables to use in this case are Xjjclc) = Y. ^+,crn'n.®a]^. 
An example is the interaction [(n) = rL^(3/10)’^. The action of the map F+ for 
this interaction is illustrated figure 1. Plotted are the pairs {4)(+Vc7), 4)(+V+V 
cr)}. The plot does not represent a function, as there are points that have several 
images. A experimental point of view can be taken, and the plot considered 
as the chaotic attractor for an undetermined dynamical system. It should then 
be analyzed to determine the dimensionality of the dynamical system and how 
many dimensions are required to embedded it smoothly(what is D for the vector 
of state observables). Only the values of the interaction that actually occur as 
a value of a configuration are plotted, and in between the points one is free to 
connect them as required for the smoothness of the map F+. 

The added smoothness and trace class of the £ operator translates into faster 
convergence towards the thermodynamic limit. If the reconstructed dynamics is 
analytic, then the convergence towards the thermodynamic limit will be faster 
than exponential [19, 20]. I will illustrate this with the polynomial-exponential 
interactions (9) with y = 1; the convergence is faster than exponential if y > 1. 
The convergence is illustrated in figure 2 for the interaction which 

requires D = 3. Plotted in the graph to illustrate the transfer matrix conver¬ 
gence are the number of decimal digits that remain unchanged as the range 
of the interaction is increased. Also in the graph are the number of decimal 
digits that remain unchanged as the largest power of iiC^ considered. The plot 
is effectively a logarithmic plot and straight lines indicate exponentially fast 
convergence. The expansion for the C operator has curvature which indicates 
that the convergence is faster than exponential. By fitting, one can verify that 
the free energy is converging to its limiting value as exp(—). Cvitanovic 
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Figure 1: The spin adding map F+ for the potential J(n,) = The action 

of the map takes the value of the interaction energy between + and the semi-infinite 
configuration [oi, oi, 0 * 3 ,. ..} and returns the interaction energy between + and the 
configuration [+, ffi, 02 , 0*3 ,.. .}. 


[19] has estimated that the Fredholm determinant of a map on a D dimensional 
space should converge as exp(—+ which is confirmed by these numerical 
simulations. 

The geometrization of spin systems strengthens the connection between sta¬ 
tistical mechanics and dynamical systems. It also further establishes the value 
of the Fredholm determinant of the £ operator as a practical computational tool 
with applications to chaotic dynamics, spin systems, and semiclassical mechan¬ 
ics. The example above emphasizes the high accuracy that can be obtained: by 
computing the shortest 14 periodic orbits of period 5 or less it is possible to 
obtain three digit accuracy for the free energy. For the same accuracy with a 
transfer matrix one has to consider a 256 x 256 matrix. This make the method 
of cycle expansions feasible for analytic calculations, and it also shows that a 
few configurations are dominating the behavior of the system. 

There are many possibilities for developing the results presented. The cycle 
expansion can be factorized [3], the approach of Feigenbaum to the renormaliza¬ 
tion group can be used [ 21 ], and phase transitions can be explored as limiting 
case where the coefficients qi^ go to 1. The results also show that there is room 
for improvement in the convergence of Monte Carlo simulations. Rather than 
carrying out the averages in the spin systems, one can study averages of the 
dynamical system where convergence is faster. 

I would like to acknowledge the hospitality of the Niels Bohr Institute and the 
Deutsche Bundesbahn where this work was carried out under the NATO/NSF 
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Figure 2: Number of digits for the Fredholm method (•) and the transfer function 
method (x). The size refers to the largest cycle considered in the Fredholm expansions, 
and the truncation length in the case of the transfer matrix. 
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